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Abstract 

Motivated by the application to spacetimes of general relativity we 
investigate the geometry and regularity of Lorentzian manifolds under cer- 
tain curvature and volume bounds. We establish several injectivity radius 
estimates at a point or on the past null cone of a point. Our estimates are en- 
tirely local and geometric, and are formulated via a reference Riemannian 
metric that we canonically associate with a given observer (p, T) -where p 
is a point of the manifold and T is a future-oriented time-like unit vector 
prescribed at p. The proofs are based on a generalization of arguments 
from Riemannian geometry. We first establish estimates on the reference 
Riemannian metric, and then express them in term of the Lorentzian met- 
ric. In the context of general relativity, our estimates should be useful to 
investigate the regularity of spacetimes satisfying Einstein field equations. 

1 Introduction 

Aims of this paper 

The regularity and compactness of Riemannian manifolds under a priori bounds 
on geometric quantities such as curvature, volume, or diameter represent 
important issues in Riemannian geometry. In particular, the derivation of 
lower bounds on the injectivity radius of a Riemannian manifold, and the 
construction of local coordinate charts in which the metric has optimal reg- 
ularity are now well-understood. Moreover, Cheeger-Gromov's theory pro- 
vides geometric conditions for the strong compactness of sequences of mani- 
folds and has become a central tool in Riemannian geometry. See for instance 
(Hll|5HZH8H9Kl6H2TH221. 
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Our objective in the present paper is to present some extension of these clas- 
sical techniques and results to Lorentzian manifolds. Recall that a Lorentzian 
metric is not positive definite, but has signature (-, +,...,+). Motivated by re- 
cent work by Anderson [2] and Klainerman and Rodnianski 1 19 1, we derive here 
several injectivity radius estimates for Lorentzian manifolds satisfying certain 
curvature and volume bounds. That is, we provide sharp lower bounds on the 
size of the geodesic ball around one point within which the exponential map is a 
global diffeomorphism and, therefore, we obtain sharp control of the manifold 
geometry. Our proofs rely on arguments that are known to be flexible and effi- 
cient in Riemannian geometry, and are here extended to the Lorentzian setting; 
we analyze the properties of Jacobi fields and rely on volume comparison and 
homotopy arguments. 

In our presentation (see for instance our main result in Theorem 1 1.1 1 at the 
end of this introduction) we emphasize the importance of having assumptions 
and estimates that are stated locally and geometrically, and avoid direct use of 
coordinates. When necessary, coordinates should be constructed a posteriori, 
once uniform bounds on the injectivity radius have been established. 

Our motivation comes from general relativity, where one of the most chal- 
lenging problems is the formation and the structure of singularities in solutions 
to the Einstein field equations. Relating curvature and volume bounds to the 
regularity of the manifold, as we do in this paper, is necessary before tackling 
an investigation of the geometric properties of singular spacetimes satisfying 
Einstein equations. (See, for instance, 121 El for some background on this sub- 
ject.) 

Two preliminary observations should be made. First, since the Lorentzian 
norm of a non-zero tensor may vanish it is clear that only limited information 
would be gained from an assumption on the Lorentzian norm of the curvature 
tensor. This justifies that we endow the Lorentzian manifold with a "reference" 
Riemannian metric (denoted by gj below); this metric is defined at a point p 
once we prescribe a future-oriented time-like unit vector T. We refer to the pair 
(p, T) as an observer located at the point p. This reference vector is necessary 
in order to define appropriate notions of conjugate and injectivity radii. (See 
Section|2]below for details.) 

Secondly, we rely here on the elementary but essential observation that, in 
the flat Riemannian and Lorentzian spaces, geodesies (are straight lines and 
therefore) coincide. Under our assumptions, we will see that geodesies associ- 
ated with the given Lorentzian metric are comparable to geodesies associated 
with the reference Riemannian metric. On the other hand, the curvature bound 
assumed on the Lorentzian metric implies, in general, no information on the 
curvature of the reference metric. As we show below, one of the main is- 
sues is to guarantee the regularity of a foliation of the manifold by spacelike 
hypersurfaces. 
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Earlier work 



Let us briefly review some classical results from Riemannian geometry. Let 
(M, g) be a differentiable n-manif old (possibly with boundary) endowed with 
a Riemannian metric g. (Throughout the present paper, the manifolds and 
metrics are always assumed to be smooth.) Denote by CB(p, r) the corresponding 
geodesic ball centered at p e M and with radius r > 0. Suppose that at some 
point p e M the unit ball 23(p, 1) is compactly contained in M and that the 
Riemann curvature bound and the lower volume bound, 

l|Rm,|| L - ( Boa)) < K Vol g CB(p, 1)) > v , (1.1) 

hold for some constants K, Vq > 0. (We use the standard notation L'", 1 < 
m < oo, for the spaces of Lebesgue measurable functions.) Then, according to 
Cheeger, Gromov, and Taylor [ 10 [ the injectivity radius Inj (M, p) at the point p 
is bounded below by a positive constant i\ = i\ (K, Vq, n), 

ln) g (M,p)>h. (1.2) 

It should be noticed that this is a local statement; for earlier work on the 
injectivity radius see ll5l [TTl[T5l . 

Moreover, Jost and Karcher l!T6l relied on the regularity theory for elliptic 
operators and established the existence of coordinates in which the metric has 
optimal regularity and are defined in a ball with radius i% = i2(K,Vo,n). Precisely, 
given e > and < y < 1 there exist a positive constant C(e, y) (depending also 
upon (K, vq, n) and a system of harmonic coordinates defined in the geodesic 
ball 23(p, z'2) in which the metric g is close to the Euclidian metric in these 
coordinates and has optimal regularity, in the following sense: 

e- £ g E < g < e E g E , 

r 1+r \\dg\\ c{ ^ r)) <C( £ ,y), re(0,i 2 ]. 

Here, C° and C are the spaces of continuous and Holder continuous functions, 
respectively. Harmonic coordinates are optimal |12| in the sense that if the 
metric is of class C ,c ''' in certain coordinates then it has at least the same regularity 
in harmonic coordinates. 

The above results were later generalized by Anderson [1J and Petersen [22 J 
who replaced the L°° curvature bound by an L"' curvature bound with m > n/2. 
For instance, one can take m = 2 in dimension n = 3 in the application to 
general relativity (since time-slices of Lorentzian 4-manifolds are Riemannian 
3-manifolds). 

It is only more recently that the same questions were tackled for Lorentzian 
(n + l)-manifolds (M,g). Anderson [2, 3] studied the long-time evolution of 
solutions to the Einstein field equations and formulate several conjecture. In 
particular, assuming the Riemann curvature bound in some domain Q 

l|Rm y llL-<p) < K (1.4) 
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and other regularity conditions, he investigated the existence of coordinates 
that are harmonic in each spacelike slice of a time foliation of M. This work by 
Anderson motivated our investigation in the present paper. 

On the other hand, motivated by applications to general relativity and 
nonlinear wave equations using harmonic analysis tools, Klainerman and Rod- 
nianski [19] considered asymptotically flat spacetimes endowed with a time 
foliation and satisfying the L 2 curvature bound 

HRm*lb<E) < K (1.5) 

for every spacelike hypersurface E. To control the injectivity radius of past null 
cones, they relied on their earlier work fTTllTBl on the conjugate radius of null 
cones in terms of Bell-Robinson's energy and energy flux, and derived in 1 19 1 a 
new estimate for the null cut locus radius. We refer to these papers for further 
details and references on the Einstein equations. Section[6]of the present paper 
is a prolongation of the work (19| . 



Outline of this paper 

The present paper establishes four estimates on the radius of injectivity of 
Lorentzian manifolds, which hold either in a neighborhood of a point or on the 
past null cone at a point. Our assumptions are formulated within a geodesic 
ball (or within a null cone) and possibly apply in a ball with arbitrary size as 
long as our curvature and volume assumptions hold. All assumptions and 
statements are local and geometric. 

An outline of the paper is as follows. In Section |2j we begin with basic 
material from Lorentzian geometry and we introduce the notions of reference 
metric and exponential map for Lorentzian manifolds. In Section|3l we state our 
first estimate (Theorem 13.11 below) for a class of manifolds that have bounded 
curvature and admit a time foliation by slices with bounded extrinsic curvature. 
In Section SJ we provide a proof of this first estimate and we introduce a 
technique that will be used (in variants) throughout this paper; we combine 
two main ingredients : sharp estimates for Jacobi fields along geodesies, and 
an homotopy argument based on contracting a possible loop to two linear 
segments. In Section[5j our second main result (Theorem l5.lt shows, under the 
same assumptions, the existence of convex functions (distance functions) and 
convex neighborhoods; this result leads us to a lower bound of the convexity 
radius. 

In Section [6j our third estimate (Theorem 16. Il l covers the case of null cones 
under the assumption that the manifold has L 2 bounded curvature on every 
spacelike slice; this provides a generalization and an alternative proof to the 
result by Klainerman and Rodnianski in 1 19 ] . 

Next, in Section we establish our principal and fourth result (stated in 
Theorem l 1 . 1 I below) which provides an injectivity radius bound under the mild 
assumption that the exponential map exp is defined in some ball and the 
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curvature Rm is bounded. Most importantly, this is a general result that does 
not require a time foliation of the manifold but solely a single reference (future- 
oriented time-like unit) vector T at the base point p. This is very natural in the 
context of general relativity and (p, T) is interpreted as an observer at the point 

V- 

Given an observer (p, T), we can define the ball Bj(Q, r) C T p M with radius 
r, determined by the reference Riemannian inner product at p, and we can 
also define the geodesic ball CBr(p, r) '■— exp p (Br(0, r)). In turn, the radius of 
injectivity Inj^(M, p, T) is defined as the largest radius r such that the exponential 
map is a diffeomorphism from £>t(0, r) onto Urip, r). Let us then consider an 
arbitrary geodesic y = y(s) initiating at p and let us ^-parallel transport the 
vector T along this geodesic, defining therefore a vector field T y along this 
geodesic, only. At every point of y we introduce the reference metric gj and 
compute the curvature norm |RmJr • This allows us to express the curvature 
bound. For the convenience of the reader we state here our main result and 
refer to Section[7]for further details. 

Theorem 1.1 (Injectivity radius of Lorentzian manifolds). Let M be a time- 
orientable Lorentzian, differentiable (n + l)-manifold. Consider an observer (p,T) 
consisting of a point p e M and a reference (future-oriented time-like unit) vector 
T e TpM. Assume that the exponential map exp p is defined in the ball £>r(0, ro) c T p M 
and the Riemann curvature satisfies 

1 

sup \Rm g \ Ty < -j, (1.6) 

y r t> 

where the supremum is over the domain of definition ofy and over every g-geodesic y 
initiating from a vector in the Riemannian ball Bt(0, ro) c T p M. Then, there exists a 
constant c(n) depending only on the dimension n such that 

ln)(M,p,T) ValJ'BT(j>,c(n)r j) 

—5- >c(n) ^ '-. (1.7) 

This result should be compared with the injectivity radius estimate estab- 
lished by Cheeger, Gromov, and Taylor [10 [ in Riemannian geometry. Observe 
that the curvature assumption | |1.6t can always be satisfied by suitably rescal- 
ing the metric tensor. It would be interesting to replace the volume term in the 

right-hand side of lO by Vol ? (S(p, r )). 

Finally, in the last two sections of this paper, we establish a volume compari- 
son theorem for future cones and generalize our main theorem to the volume of 
a future cone (Section|8]l, and we briefly discuss the regularity of the Lorentzian 
metric in harmonic-like coordinates, and present a generalization to pseudo- 
riemannian manifolds (Section|9]l. 
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2 Preliminaries on Lorentzian geometry 



Basic definitions 

It is useful to discuss first some basic definitions from Lorentzian geometry, 
for which we can refer to the textbook by Penrose l20l . Throughout this pa- 
per, (M,g) is a connected and differentiable (n + l)-manifold, endowed with a 
Lorentzian metric tensor g with signature (— , +,...,+). To emphasize the role 
of the metric g or the point p we use any of the following notations 

g p (X, Y) = (X, Y) gf = <X, Y) g = (X, Y) p 

for the inner product of two vectors X, Y at a point p e M; we sometimes also 
write |X|| instead of g p (X, X). Recall that the tangent vectors X e T p M are called 
time-like, null, or spacelike depending whether the norm g p (X, X) is negative, 
zero, or positive, respectively. Vectors that are time-like or null are called causal. 

The time-like vectors form a cone with two connected components. The 
manifold (M, g) is said to be time-orientable if we can select in a continuous 
way a half-cone of time-like vectors at every point p. The choice of a specific 
orientation allows us to decompose the cone of time-like vectors into future- 
oriented and past-oriented ones. The set of all future-oriented, time-like vectors 
at p and the corresponding bundle on M are denoted by T p M and T + M, respec- 
tively. We also introduce the bundle T+M consisting of elements of T + M with 
unit length. 

By definition, a trip is a continuous curve y : {a, V) M made of finitely 
many future-oriented, time-like geodesies. We write p « q if there exists a trip 
from p to q. A causal trip is defined similarly except that the geodesies may be 
causal instead of time-like, and we write p < q if there exists a causal trip from 
p to q. 

The set 3 + {p) := [q e M / p « is called the chronological future of the 

point p, and 3~(p) := uj £ M / q « pj is called the chronological past. The 
causal future and past are defined similarly by replacing << by <. The future 
or past sets of a set S c M are defined by 

3^):= \jF(p), r(S):=|Jr(p), 

peS peS 

and one easily checks that ± (S) are open, but that / ± (S) need not be closed in 
general. 

A future set F c M by definition has the form F = 3 + (S) for some set S c M. 
Similarly, a past set satisfies F = J~(S) for some S. A set is called achronal if no 
two points are connected by a time-like trip. Observe that a set can be spacelike 
at every point without being achronal and that an achronal set can be null at 
some (or even at every) point. A set B cz Mis called an achronal boundary if it 
is the boundary of a future set, that is B = d3 + (S) = 3 + (S) \ 3 + (S) for some ScM. 
One can check that given a non-empty achronal boundary the manifold can be 
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partitioned as M = P U B U F, where B is the boundary of both F and P and, 
moreover, any trip from p e P to q e F meets B at a unique point. Observe also 
that any achronal boundary is a Lipschitz continuous n-manif old. For instance, 
in Section|6]below, we will be interested in the geometry of past null cones, that 
is the sets d3~(p) for p e M. 

Given an arbitrary achronal and closed set S C M, we define the (future or 
past) domains of dependence of S in M by 

D ± (S) := [peM/ every future (resp. past) endless trip containing p meets S , 
D(S) := D"(S)UB + (S). 

Observe that domains of dependence are closed sets. Next, define the (future 
or past) Cauchy horizons 

!K*(S) := [p e D ± (S) / D*(p) n © ± (S) = 0} = D ± (S) \ 3 T (D ± (S)), 
JC(S) := 5T(S) U IK + (S). 

For instance, the future Cauchy horizon is the future boundary of the future 
domain of dependence of S. One can check that the Cauchy horizons are closed 
and achronal sets, with <?£> + (S) = W + (S) U S and dD(S) = 5f(S). 

Finally, a (future) Cauchy hypersurface for M is defined as an achronal (but 
not necessarily closed) set S satisfying D + (S) = M. For instance, it is sufficient 
for S to be smooth, achronal, spacelike and such that every endless null geodesic 
meet M. 

Reference metric 

As explained in the introduction one should not use the Lorentzian metric to 
compute the norm of a tensor since the Lorentzian norm may vanish even 
when the tensor does not. This motivates the introduction of a "reference" 
Riemannian metric associated with a time-like vector field, as follows. 

Let T be a future-oriented, time-like, unit vector field, satisfying therefore 
g p (T, T) = -1 at every point p. We refer to T as the reference vector field pre- 
scribed on M. Introduce a moving frame E a (a = 0, 1, . . . , n) defined in M, that 
is, E a is an orthonormal basis of the tangent space at every point and consists 
of the vector eo = T supplemented with n spacelike unit vectors e; (j = 1, . . . , n). 
Denoting by E a the corresponding dual frame, the Lorentzian metric tensor 
takes the form 

g = i h f S E a ®Ef 1 , 

where r\ a p is the Minkowski "metric". This decomposition suggests to consider 
the Riemannian version obtained by switching the minus sign in rjoo = _ 1 into 
a plus sign, that is 

gT :=b a pE a ®Ef', 

where 5 a B is the Euclidian "metric". Clearly, gj is a positive definite metric; it 
is referred to as the reference Riemannian metric associated with the frame E a . 
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For every p 6 M, since T p is time-like, the restriction of the metric g p to 
the orthogonal complement {Tp} 1 - c TpM is positive definite, and the reference 
metric can be computed as follows: if V = a T p + V and W = bT p + W with 
V, W e {Tp}- 1 , then 

g T , p (V,W)=ab + g p (V',W). 
In the following, we use the notation 

g T/V {V, W) = (V, W)t, p , gi, v {V, V) = \V\ 2 Zp 

for vectors; the norm of tensors is defined and denoted similarly. 

In contrast with the Lorentzian norm, the Riemannian norm |A|r, p of a tensor 
A at a point p e M vanishes if and only if the tensor is zero at p. Moreover, 
as long as T remains in a compact subset of the bundle of half-cone T + M, the 
norms associated with different reference vectors are equivalent. 

The reference Riemannian metric also allows one to define the functional 
norms for Lebesgue and Sobolev spaces of tensors defined on M (as well as 
on submanifolds of M), allowing us for instance to view L 2 (M, gj) as a Banach 
space. In particular, we will use later the L 2 norm of a tensor field T on M 
restricted to an hypersurface E: 

\Nh\\V(Z, g r) := §m\dV ZigT , 

where dV^g? is the volume form induced on E by the reference Riemannian 
metric. The functional norm above depends upon the choice of the vector field 
T, but another choice of T would give rise to an equivalent norm (provided T 
remains in a fixed compact subset). Observe in passing that the volume forms 
associated with the metrics g and gj coincide, so that the spacetime integrals 
of functions in (M, g) or (M, gj) coincide; for instance, the volume Vol ? (A) and 
Volg T (A) of a set A c M coincide. 

Furthermore, we observe that in order to define the reference metric gj at a 
given point p, it suffices to prescribe a future-oriented time-like unit vector T at 
that point p only; it is not necessary to prescribe a vector field. In the situation 
where the reference metric need only be defined at the base point p, we refer to 
T as the reference vector (rather than vector field) and we refer to (p, T) e T+M 
as the observer at the point p. This will be the standpoint adopted for our main 
result in Section[7]below. 

Exponential map 

On a complete Riemannian manifold the exponential map exp p : T p M — » M at 
some point p e Mis defined on the whole tangent space T p M and is smooth. 
For sufficiently small radius r the restriction of exp p to the ball B ?ji (0, r) c T p M 
(determined by the metric g p at the point p) is a diffeomorphism on its image. 
The radius of injectivity at the point p is defined as the largest value of r such 
that the restriction exp „\B gp (o,r) is a global diffeomorphism. 
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In the Lorentzian case, the exponential map is defined similarly but some 
care is needed in defining the notion of radius of injectivity. First of all, if 
the manifold is not geodesically complete (which is a rather generic situation, 
as illustrated by Penrose and Hawking's incompleteness theorems [14 1), the 
map exp p need not be defined on the whole tangent space T p M but only on 
a neighborhood of the origin in T p M. More importantly, the Lorentzian norm 
of a non-zero vector may well vanish; consequently, the radius of injectivity 
should not be defined directly from the Lorentzian metric g. The definition 
below depends on the prescribed Riemannian metric gT, f , at the point p. 

Definition 2.1. The conjugate radius Conj^(M, p, T) of an observer (p, T) e T+M is 

the largest radius r such that the exponential map exp p is a local diffeomorphism from 
the Riemannian ball Bt(0, r) = B gTi> (0, r) c T p M to a neighborhood ofp in the manifold 
M. Similarly, the injectivity radius Inj „(M,p, T) of an observer (p, X) 6 T+M is the 
largest radius r ssuch that the exponential map is a global diffeomorphism at every 
point of the ball B T (0, r). 

When a vector field T is prescribed on the manifold (rather than a vector at 
the point p), we use the notation Inj (M,p, T) instead of Inj ? (M, p, T p ). The radii 
Conj (M,p, T) and Inj (M,p, T) are essentially independent of the choice of the 
reference vector, as long as it remains in a fixed compact subset of T p M. 

We also need the notion of injectivity radius for null cones. Given a point 
p 6 M and a reference vector T e T p M, we consider the past null cone at p, 

N- :=\XeT p Mj g p (X, X) = 0, g p (T, X) > oj, 

which is defined a subset of the tangent space at p. Denote by 

B^(0 / r) = B^(0 / r):=B gT ^0,r)nN- 

the intersection of the Riemannian gr p -ball with radius r and the past null cone, 
and by 

K"(p) := dT{p) 

the past null cone at p. 

Consider now the restriction of exp p to the past null cone, denoted by 

exp^ : B£(0, r) c N~ -» JT(p) c M, 
which we refer to as the null exponential map. 

Definition 2.2. The (past) null conjugate radius Null Conj (M, p, T) of an observer 
(p, T) e T^M is the largest radius r such that the null exponential map exp^ is 
a local diffeomorphism from the punctured Riemannian ball Bj(0, r)\{0] c T p M to a 
neighborhood ofp in the past null cone. The null injectivity radius Null Inj (M, p, T) 
of an observer (p, T) 6 T^M is defined similarly by requiring the map exp^ to be a 
global diffeomorphism. 
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3 Lorentzian manifold endowed with a reference 
vector field 



A first injectivity radius estimate 

From now on, we fix a reference vector field T which allows us to define the 
Riemannian metric gj and compute the norms of tensors. We begin with a set 
of assumptions encompassing a large class of Lorentzian manifolds with L°° 
bounded curvature and we state our first injectivity estimate, in Theorem 13. II 
below. The forthcoming sections will be devoted to further generalizations and 
variants of this result. 

We fix a point p e M and assume that a domain Q c M containing p 
is foliated by spacelike hypersurfaces E< with future-oriented time-like unit 
normal T, 

Q= |J E t . (3.1) 

(€[-1,1] 

The positive coefficient n is defined by the relation 4 = n T, or 




In the context of general relativity, n is the proper time of an observer mov- 
ing orthogonally to the hypersurfaces, and is called the lapse function. The 
geometry of the foliation is determined by this function n together with the 
Lie derivative &rg- The latter is nothing but the second fundamental form, or 
extrinsic curvature, of the slices Ef embedded in the manifold M. 

We always assume that the geodesic ball "Bzgip, 1) c Eo (determined by the 
induced metric g|z ) is compactly contained in Eo. We introduce the following 
assumptions: 

e~ K ° <n<e K ° in Q, (Al) 
\H T g\ T < Ki in Q, (A2) 
IRrrtjIr < K 2 in Q, (A3) 

Vol gko CB Eo (p,l))>»o, (A4) 

where K.q,K\,K2 and Vq are positive constants. Observe that Assumption (A4) 
is a condition on the initial slice only; together with the other assumptions it 
actually implies a lower volume bound for every slice of the foliation. 
We will prove: 

Theorem 3.1 (Injectivity radius of foliated manifolds). Let Mbe a differentiable 
manifold endowed with a Lorentzian metric g satisfying the regularity assumptions 
(A1)-(A4) at some point p and for some foliation (|3.11 . Then, there exists a positive 
constant z'o depending only upon the foliation bounds Kq,K\, the curvature bound K2, 
the volume bound Vo, and the dimension of the manifold such that the injectivity radius 
at p satisfies 

Inj (M,p,T)>z . 
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The following section is devoted to the proof of this theorem. Observe 
that the conditions (A1)-(A4) are local about one point of the manifold and are 
stated in purely geometric terms, requiring no particular choice of coordinates. 
Of course, the conclusion of Theorem |3 . 1 I hold globally in M if the assumptions 
(A1)-(A4) hold also globally at every point of the manifold. Our assumptions do 
depend on the choice of the time-like vector field T, but the dependence of the 
constants arising in (A1)-(A4) should not be essential; however, it is conceivable 
that, when applying this theorem in a specific situation a quantitatively sharper 
estimate would be obtained with a choice of an "almost Killing" field, that is 
a field T corresponding to a "small" Lie derivative Urg. Later in Section a 
more general approach is presented in which the vector field T is constructed 
from a single vector prescribed at the point p. 

Basic estimates on the reference metric 

To establish Theorem |3.1| it is convenient to introduce coordinates on Q, chosen 
as follows. Fix arbitrarily some coordinates (x l ) on the initial slice Eo- Then, 
transport these coordinates to the whole of Q along the integral curves of the 
vector field T. This construction generates coordinates (x a ) on Q such that x° = t 
and the vector d/dt is orthogonal to d/dx', so that the Lorentzian metric takes 
the form 

g = -n 2 dt 2 + gijdtfdxK (3.2) 

where n is the lapse function and gjj is the Riemannian metric induced on the 
slices Ef . The reference Riemannian metric in the domain Q then takes the form 

g T = n 2 dt 2 + gijdx'dx', (3.3) 

and the Riemannian norm of a vector X has the explicit form: gj(X,X) := 
n 2 X°X° + X'Xj. 

We want to control the discrepancy between the reference Riemannian met- 
ric gj and the original Lorentzian metric g, as measured in the connections 
V and V ?T and the curvature tensors Rm and Rm Sr Clearly, these estimates 
should involve the constants arising in (A1)-(A4). Consider the general class 
of metrics 

g:= fdf + gijdx'dx', (3.4) 

which allows us to recover both the Lorentzian (/ = -n 2 ) and the Riemannian 
(f = n 2 ) metrics. 

In view of the expressions of the Christoffel symbols and the Riemann 
curvature _ 

fy . = l~ys _ dg a p \ 

«P 4 2 g V dx« dx?> dx b >' 

dl\, dl c 

af}5 ■ faa J x p m l (id fir] ab' 
Rafiyd ■= gyi^a^ R afi ■= Rayfib^ 
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we compute explicitly the Christoffel symbols associated with the metric g, 

(3-5) 



~ \_dj_ f 0= J_d£ y 0= _L d JH 

00 2f dt' Qi 2fdx 1 ' { i 2/ dt 



00 ~ 2 g dx 1 ' i0 ~ 2 s dt' l '~ v' 



as well as the non-trivial curvature terms 



„ „ 1 (dg^dgjl dgudgjk, 

Riju-Rijki- Tf {-sf^f--df^fh 



rp JH.A. (l^ d J±\ + V-**i*, - r p ( V^) 

V # toV2 g dt > 2 g dt gqk ' l > kK 2 g dt' 

V4/ x At*' <?f 2* * 2/ ^x'' 
~ l(„fd \ „(d \\ 1 /df dgn df dgijs 

100 dA 2 g dxii dt^2 g dt > 'A 2 g dxi> 

dt dt i 2f dt^ dt > 2 s dxt2fdx 



and 



r - 1 Cv v f + ^ U 1 «*? ^ + 1 ^ ^ + 1 d f d f 

Ri °i° - "2i w + -dF) + ^r^r + i/^-^r + 17^*7- 

By applying the formulas above to both metrics g,gr we estimate the 
Christoffel symbols, as follows. Recall that the difference T y a „ — T r ^ ~ can be 
regarded as a tensor field on M, so that the following (Riemannian) norm 
squared is a scalar field on the manifold M: 

iv*. - v| 2 T := |r, T - n 2 T = cr^ - 1-) (T^ r - rf Y )g T , aa , gf 

We need also the expression of the Lie derivative of g along the vector field T. 
By a direct computation from l|3.2l l we obtain 

(£rg)oo = 0, (Zrgh = \% ^Tg)ij = ^. (3.6) 

n ox' n dt 

Lemma 3.2 (Levi-Cevita connection of the reference metric). Suppose that g 
satisfies Assumptions (A1)-(A2). Then, the covariant derivative of the Lorentzian and 
Riemannian metrics are comparable, precisely 

|V fr - V| T = n 2 \L T g\ 2 < e 2K "Kl =: K 3 . 
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Proof. In view of <|3 .5b the difference T gT - T depends essentially upon the terms 

jg and ^ which precisely appear in the expression of the Lie derivative j3.6\ . 
We omit the details. □ 

It is important to observe that the difference between the curvature tensors 
can not be similarly estimated, and that this is one of the main difficulties to 
deal with in the present work. 

For future reference we provide here the expressions of certain curvature 
coefficients of g and gj in terms of (first-order derivatives of) the lapse function 
n and the induced metric g»; 

l(„,d v ( d ,\ 1 ( dn 2 dgj, dn 2 dgij, 

d _ 1 /'Y7Y7/2x dZ Sii\ 1 V n d Sip d Sn 1 dfl 2 dgij 1 dtl 2 dn Z 

and 



KTAjB-Hijti An2 { ft ft dt dt )' 

l /v7V7 . 2 , ^ 2 gyx 1 dgtpdgjc, 1 ^n 2 ^; 1 5« 2 5n 2 

where R z .,, denotes the induced curvature tensor on the time slices E = E f . 

4 Derivation of the first injectivity radius estimate 

In this section we provide a proof of Theorem l3.ll 

Step 1. Radius of definition of the exponential map. First of all, we note 
that the injectivity radius of the Riemannian metric g\z, induced on the initial 
hypersurface Lq = t _1 (0) is controled, as follows. Using Assumptions (A3) and 
(A4), we see that the Riemann curvature of the metric g\z a is bounded and the 
volume of the unit geodesic ball Vol ? | E(j (23e (P/ 1)) is bounded below. Therefore, 
according to IflOl , there exists a constant i\ - h(K2,vo) such that the injectivity 
radius of g\z at the point p is i\ at least: 

Inj^Ecp) ^ h- 

Moreover, according to [16J we can also assume that i\ is sufficiently small so 
that, given any e > there exists a coordinates (x a ) defined in a ball with definite 
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size near p, with x a (p) = 0, such that the metric g\z is close to the n-dimensional 
Euclidian metric = 5u (in these coordinates). More precisely, on the initial 
slice Eo we have 

e~ e 5ij < gij (0, x 1 ,..., x") < fy, (x\ x") e B E ,(0, k), 

where we have set B F (0, r) := {(x 1 ) 2 + . . . + (x") 2 < r 2 } c W. The latter can be 
regarded as a subset of Eo by identifying a point with its coordinates, and we 
also use the notation 23 p (p, r) for this Euclidian ball. 

We can next introduce some coordinates (x a ) = (t,x ] ) on the manifold, by 
propagating the coordinates (x 1 ) chosen on Eo along the integral curve of the 
vector field T. This construction allows us to cover the domain Q. From 
Assumption (A2) (together with (Al) and i3.6l l) we deduce that the induced 
metric on each slice of the foliation is comparable with the n-dimensional 
Euclidian metric in some time interval \—iz, i 2 ], that is 

(e- E -Ki 2 ) 5ij < glj (x) < (e c +Ki 2 ) b ip 
x = (t,x 1 ,...,x n )e H2,;'2]xB P (0,/i), 

for some K > depending only on Kq, K\, K 2 . 

We then restrict attention to a smaller radius i 2 = h(Ko, K\, K 2 ) < i\ such that 
e~ f - Ki 2 > 0, and we pick up c\ > sufficiently large so that e~ Cl < e~ £ - Ki 2 < 
e e + Ki 2 < e Cl . In turn, in view of Assumption (Al) on the lapse function n 
and of the expression (|3.3b of the reference Riemannian metric gj, the above 
inequalities imply that the reference Riemannian metric gj is comparable to the 
(n + l)-dimensional Euclidean metric: 

e _C2 Sop < gT, a p < e C2 b ap x = (t,x\...,x n )e [-i 2 , i 2 ] x B E , (0, i 2 ) 

for some constant c 2 > C\ depending upon C\ and Kg. 

Introducing on the manifold the (n + l)-dimensional Euclidian metric E 
(which we define in the constructed coordinates (x a ) and is, of course, indepen- 
dent of the point on the manifold) and the corresponding Euclidian metric ball 
23e(P/ h), we have established 

e~ Cl g E < gT, q < e Cl gE, q 6 S E (p, i 2 ). (4.1) 

In the following we use the notation |X|e for the Euclidian norm of a vector X. 

Our first task is to determine the radius of a ball on which the exponential 
map is well-defined. This radius depends upon the reference vector field T. 
Let y : [0, So] — > M be a geodesic associated with the Lorentzian metric g and 
satisfying y(0) = p. Assume that this geodesic is included in the Euclidian ball 
^E(P/h) ( m which we have a good control of the metric gj)- Obviously, we 
have 

(y'(s),y'(s)) g = </(0),y'(C%, s e [0,s ]. 
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On the other hand, to determine the length of y'(s) with respect to the reference 
metric gj, we proceed as follows: 

\l(y'(s),y'(s)) T \ = |v r , /(s) (g r (y'(s),y'(s)))| = 2 |<V ?T//(s) /(s),/( S )) r | 

= 2|<(V ?T -V ? V (s) /(s),y'(s)) T | 
<2|V, T -V ? | T |/(s)| 3 T . 

So, by Lemma l3~2l 1 (s)|^,| < 2X3 |y'(s)|^, and, in consequence, 

By integration of the above inequality and provided s is small enough so 
that 2s K 3 |y'(0)| r < 1, we see that 

i|/(0)lr<lr'(s)lT<2|/(0)| T . (4.2) 
In view of (|4.1| | this implies 

^ 1/(0)12 < |/(s)| E < 2 ^ |y'(0)| E . (4.3) 

These inequalities hold for all s e [0, 1/(2K 3 |y'(0)| T )] as long as y(s) e B E (p, i 2 ). 
In particular, by restricting attention to geodesies whose initial vector has unit 
Euclidian length, |y'(0)| E = 1, we see that y([0, r 2 ]) c T> E (p, i 2 ) where r 2 
i 2 e~ Cl 12. In turn, this establishes that the exponential map at the point p is 
well-defined on the ball B E (0, r 2 ) with a range included in the geodesic ball 
£ E (p,z' 2 ). 

Step 2. Conjugate radius estimate. Our second task is to determine a ball on 
which the exponential map is a local diffeomorphism, and we therefore need 
to control the length of a Jacobi field along a geodesic. Let y : [0, r 2 ] — » M be a 
^-geodesic satisfying y(0) = p and |y'(0)|e = 1. By the discussion in Step 1 we 
already know that the curve y lies in IS^p, i 2 ) and that max se [o /r j |y'(s)|r < 2 e 2ci . 
Given an arbitrary Jacobi field along y, } = }(s), satisfying 

r(s) = -Rm(/(s),yW(s), 
/(0) = 0, |/'(0)| r = 1, 

we need to control its Riemannian length F(s) := |/|r(s), as stated in l|4.7|l below. 

Let [0, So] be the largest subinterval of [0, r 2 /4] in which the inequality |/|r < 1 
holds. Using the equation satisfied by the Jacobi field and taking into account 
the curvature bound (A3), we deduce that, in the interval [0, So], 

|^<V / 7,V / /> T | =2|(V gT// V / /,V / 7) T | 

< 2 |V ?r - V g \ T \y'\T NyJ\ 2 T + 2K 2 \y%\J\ T \V y ]\ T . 
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With ( 14.2b and the co variant derivative estimate in Lemma [3721 we obtain 

|^|V / /| r |<4K 3 |V / /| r + 8JC 2 . (4.4) 

We can next integrate l|4.4t over an arbitrary interval [0,s] c [0,So], use the 
initial condition on the Jacobi field, and obtain 

1 + ^ (1 - e~^ s ) < |VWIr < 1 + ^ - 1). 

Assuming that r 2 is small enough so that ^ (1 - e~ 4K3S ) < 1 /2 and ^ (e 4K;)S - 1 ) < 
we infer that 

\ < |V r /| T < 2. (4.5) 

Hence, using this inequality and Lemma [3~72l we find £\J\ T < 2 + 2X 3 < 1, so that 

F(s) = |/| r (s) < (2 + 2K 3 ) s < (2 + 2K 3 ) r 2 . (4.6) 

Further assuming that (2 + 2K 3 ) r 2 < 1 we conclude that s = r 2 . 
Next, we want to improve the rough estimate ( 14.6b . Since 

J s WylJ)T = W gT ,yV r h1)T + (V gr , Y l V r /) T 

then by substituting the previous estimates of |/|r(s) and \Vy /|r(s) and perform- 
ing similar calculations as above, we get 

e~ C3 < j s (V r J,J) T <e c i 
for some constant c 3 > 0. By integration this implies 

e~ C3 s < (V r J,J) T < e C3 s 
and we arrive at the following lower bound for the norm of the Jacobi field: 

l<V,././>r| e -«s „ 

m> --W7nr > - — > -' s 

for some C4 > 0. 

On the other hand, using again the above estimates we have 



^F<^{(V gT , r J,Jh + K 3 F 2 ) 



— F < - 
ds ~ F 

S 



< — (e C3 s + K 3 (2 + 2K 3 ) 2 s 2 ) < e C5 

for some constant C5 > 0. This leads to the upper bound 

F(s) < e Cs s. 
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In summary, we have established that the norm of the Jacobi field is com- 
parable with s: 

e~ Ci s < F(s) < e C5 s, s e [0, r 2 \. (4.7) 

By the definition of Jacobi fields these inequalities are equivalent to controling 
the differential of the exponential map, that is for s e [0, r 2 ] 

e~ Ci \W\ T < |dexp ps/(0) (W)| r < e c * \W\ T . 

We conclude that the pull back of the reference metric to the tangent space at p 
satisfies 

e~ Ci g T ,p < (exp p ) g T < e C5 g T , P 

in the ball B T (0, r 2 ) c T p M. 

In particular, since the conjugate radius of the Lorentzian metric is precisely 
defined from the reference Riemannian metric, these inequalities show that the 
conjugate radius of the exponential map is r 2 at least. 

Step 3. Injectivity radius estimate. We are now in a position to establish 
that Inj (M,p, T) > r 3 := r 2 e~ C2 /4. We argue by contradiction and assume that 
y\ : [0,si] — > M and y 2 : [0,s 2 ] -» M are two distinct g-geodesics satisfying 
max(si,s 2 ) < r 3 and 

yi(0) = 7z(0) = p, |/ x (0)|r = |/ 2 (0)| r = 1, 
yi(si) = y 2 (s 2 ) =: q. 

We will reach a contradiction and this will establish that the injectivity radius 
is greater or equal to r 3 (as can be checked by using the fact that the exponential 
map is at least a local diffeomorphism). 

By Step 1 we know that y\,y 2 c H E (p,2e 2c2 r 3 ). By concatenating these two 
curves, we construct a geodesic loop containing p, 

y = y- 1 Uy 1 : [0,s a + s 2 ] -» S E (p,2e 2c2 r 3 ), 

which need not be smooth at p or q. Since y is contained in the image of the ball 
f>r(P/ r 2 ) under the exponential map, we can define an homotopy of y with the 
origin (x = 0), by setting (in the coordinates constructed earlier) 

r e (s) = ey(s), ££[0,1]. 

The curves T t : [0,si + s 2 ] -> S E (p,2e 2c2 r 3 ) satisfy 

r £ (0) = r £ ( Sl + s 2 ) = p, r ([o,i]) = p, ii = y. 

Moreover, we have \T' £ (s)\ E < s2e 2c2 < le 202 and thus |r' £ (s)|r < 2e 3c2 . In 
particular, the gr-lengths (computed with the reference metric) of the loops T £ 
are less than 

L(r f ,g r )<2e 3c2 r 3 = |. 
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Since the exponential map is a local diffeomorphism from the ball Bj(0, rj) C 
TpM to the manifold, and in view of the estimate < t4.8b on the exponential map, 
it follows that all the loops T £ can be lifted to the ball £>t(0, r-i) in the tangent 
space with the same origin 0. Consequently, we obtain a continuous family of 
curves f £ : [0, si + s 2 ] — » T p M satisfying 

r £ (0) = 0, ee[0,l]. 

At this juncture we observe that, since T £ (si + s 2 ) (for e e [0,1]) all cover the 
same point p and since the curve To is trivial and the family is continuous, 

? £ ( Sl + S 2 ) = 0, £€[0,1]. 

It remains to consider the lift of the original geodesic loop y: under the lifting 
the geodesies 71,72 are sent to two distinct line segments (with respect to 
the vector space structure) originating at the origin which obviously do not 
intersect. This is a contradiction and we conclude that, in fact, Inj g (M,p, T) > 
as announced. This completes the proof of Theorem 13. II 

5 Convex functions and convex neighborhoods 

We establish now the existence of convex functions and convex neighborhoods 
in M. Let us recall first some basic definitions. A function u is said to be 
geodesically convex if the composition of u with any geodesic is a convex 
function (of one variable). A set Q' c Q" is said to be relatively geodesically 
convex in O" if, given any points p,q € Q' and any geodesic (segment) y from p 
to q contained in Q", one has y C Q'. A set Q' is said to be geodesically convex 
in Q" if Q' is relatively geodesically convex in Q" and, in addition, for any p, q, 
there exists a unique geodesic y connecting p and q and lying in Q' . 

We denote by dj the distance function associated with the reference Rie- 
mannian metric gj. 

Theorem 5.1 (Existence of geodesically convex functions). Let (M, g) be a differ- 
entiable (n + \)-manifold endowed with a Lorentzian metric g, satisfying the regularity 
assumptions {Al)-{AA)for some point p 6 M and some future-oriented, unit, time-like 
vector field T, and let gj be the reference Riemannian metric associated with Then, for 
any e 6 (0, 1) there exists a positive constant rg depending only upon e, the foliation 
bounds Ko, K\, the curvature bound K2, the volume bound vq, and the dimension of the 
manifold and there exists a smooth function u defined on Hjip, tq) such that 

(l-e)d T (p,-) 2 <u<(l + £ )d T (p,-) 2 , 
(2-£)g T <V 2 u<(2 + £)g T . 

Hence, the function u above is equivalent to the Riemannian distance func- 
tion from p and is geodesically convex for the Lorentzian metric. In the proof 
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given below, the function u is the Riemannian distance function associated with 
a new time-like vector field (denoted by N in the proof below). The following 
corollary is immediate and provides us with a control of the radius of convexity, 
which generalizes Whitehead theorem from Riemannian geometry 11231161. 

Corollary 5.2 (Existence of geodesically convex neighborhoods). Under the as- 
sumptions of Theorem 15.11 for any < r < r^ there exists a set Q,- c Q which is 
geodesically convex in T>T(p,2ro) and satisfies 

exp p (B T (0, r)) c Q r c exp p (B T (0, (1 + 5)r)). 

Moreover, one can always choose Q r so that 

S T (p,r)cG r cB T (p,(l + 5)r) ( 

where Ht(P, r ) is the geodesic ball determined by the reference Riemannian metric. 

Proof of Theorem 15711 Step 1. Synchronous coordinate system. Given e > 0, by 
applying the injectivity radius estimate in Theorem l3.1l to points near p, we see 
that there exists a constant rg depending on Kq, K\, K2, Vq, e, n such that for any 
q e H>t(p, 2ro) the injectivity radius at q is 2ro at least, and we can assume that 

e~ E g T , q < (exp q )*g T < e £ g TiV B T (0, r ) c T,M, q e S r (p, 2r ). 

Let y = y(s) be the (backward) time-like geodesic satisfying y(0) = p and 
y'(0) = — T p , and consider the (past) point q :- y(ro/2). The future null cone 
at q with radius ro (the orientation being determined by the vector field T) is 
defined by 

Q(ro) := [V 6 T,M / 1 < r , < 0, (V, T) > oj. 

Observe that the ^j-length of y between p and q is approximatively r /2 and that 
the norm \y'\r is almost 1, while |y'(^)l| = 1 an£ i ( — 7'/ > 0- By the injectivity 
radius estimate in Theorem 13.11 the exponential map at q is a diffeomorphism 
from Cq(ro) onto its image which, moreover, contains the original point p. 
Next, introduce the set of vectors that are "almost" parallel to T: 

Q(r , e) := [v 6 T q M.j \V\ T , q < r , (V, T> ft > 0, > 1 - e}. 

The notation c(e) > is used for constants that depend only on Ko, K\, K2, Vq, n, e 
and satisfy lim £ ^o c(e) = 0. We claim that there is constant c(e) > such that 

£ r (p,c(£)r )cexp 9 (C ? (ro,£)). (5.1) 
Actually, we have 23r(p,c(e)r ) c Uriq, (5 + c(e))ro), hence 

S T (p,c(e)r ) c exp 9 (B T (0,(i + c(e))r )). 
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Since the metrics gj,o and gj jq are comparable (under the exponential map at q) 
we see that geodesies a connecting q and points of 23r(p, c(e)ro) make an angle 
< c(e) with -y'(q) at the point q (as measured by the metric gr#)- By reducing 
the constant c(e) if necessary, the claim is proved. 

Let t be the Lorentzian distance from q: it is defined on exp i; (C ? (fo)) and is 
a smooth function on exp (C 9 (ro)) \ (p}. Using the claim d5.lt we deduce that t 
is smooth in the ball l Br(p, c(e)ro) and satisfies 

1 1 
(- - c(e)J r Q <%<(- + c(e)J r in the ball B r (p, c(e)r ). (5.2) 

It is clear also that 

|Vt|| = -1, V 2 t(Vt, ■) = 0. 

We now introduce a new foliation. Let (z 1 ) be coordinates on the level 
set hypersurface t = x(p), and by following the integral curves of the (unit, 
time-like) vector field 

N:= Vt 

let us construct coordinates (z") with Zo := t in which the Lorentzian metric g 
takes the simple form 

g = -(dz°) 2 + g,jdz'dzi. 

Let gN = (-,-)n be a (new) reference Riemannian metric based on the vector 
field N. 

By Lemma [3.2l using the equation satisfied by (future) g-geodesics a we see 
that 

|— log|a'(T)|| < X 3 r - 

(Recall that we allow ro to depend upon e.) This inequality shows that the 
vector field N makes an angle < c(e) with T, everywhere on exp (? (C £; (ro, e)). 
From this, we conclude that the two metrics are comparable: 

(1 - c(e)) g T <gN<(l + c(e)) g T in the cone exp^C^ro, e)). 



Step 2. Hessian comparison theorem and curvature bound for the reference 
metric gN- Since p e exp ? (C ? (ro)), let a : [0, x(p)] — > M be the future time-like 

geodesic connecting q to p, and let V be the Jacobi field defined along a such 
that 

V(Q) = 0, V(T(p)) = V, 
where V e T p M satisfies the orthogonality condition (Vx, V) = 0. Then we have 

-v 2 t(v, v) = -<y, v Vt v) = (v, Vj v> 



■X 



■t(p) _ _ _ 

<V a V, V i V>. - R(ff', y, a', V) =: I(V, V). 
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Recall that in the absence of conjugate points Jacobi fields minimize the index 
form I(V, V) among all vector fields with fixed boundary values. By applying a 
standard comparison technique from Riemannian geometry on the orthogonal 
space (Vt)" 1 (on which the Lorentzian metric induces a Riemaniann metric) we 
control the Hessian of t in terms of the curvature bound K2 : 

V*2(l + c(e)) VK 2 (1 + C ( £ )) 

g |(Viy-< (-V z t)| ( vt)^ < ; g I(VtK ■ (5.3) 



tan tJK 2 (1 + c(e))T tanh ^(1 + c(e))x 

Since — V?.x = we deduce from 05.3b that 

— < -rr— < m the cone exp (C q (r )). (5.4) 

Combining ||5.4| | with the curvature formulas derived in Section 3, i.e. 

"if* ~ V + 4 V dx dr dx dT )' 

d _ ld 2 gij 1 dgiphh 
Ri0 '°--2l^ + l g ' 

|^| <§ onexp 9 (C 9 (r )). (5.5) 

Finally, relying on the formulas for the curvature of the reference Rieman- 
nian metric gN, we obtain 

C 

|Ri%J N < — on exp 9 (Q(r )). 

(Observe that, as could have been expected, the upper bound blows-up as one 
approach the point q which is the base point in our definition of the distance.) 
In particular, this implies the following curvature bound near the point p: 

|Ri% N \ N < Cr~ 2 on the ball S T (p, c(e)r ). 



we conclude that 



Step 3. Constructing geodesically convex functions. Since the metrics gj arid 
gN are comparable, the volume ratio (l/rjj +1 ) Vol i r N 23]v(p,c(£)ro) is uniformly 
bounded (above and) below. By applying the theory for Riemannian metrics 
in [10[, the injectivity radius of the metric gN is bounded from below by c(e)ro. 
Let 

u(x) := d gN (p,x) 2 
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be the (square) of the distance function associated with the Riemannian metric 
gjv, which is a smooth function defined on the geodesic ball CB;v(p, c(e)ro). By 
the standard Hessian comparison theorem for Riemannian manifold we have 

pu<(2 + e) g NA p on the ball £ N (p, c(e)r ). 

In terms of the original Lorentzian metric g, the Hessian of the function u is 

Since \T gN - T\ N < C sup |^| < C by the estimate <15.4| | and since also \Vu\ N < 
2 d gN on CB N (p, r ), we conclude that 

(2 - e) g NA p > VaV^M > (2 + e) gN, a p in the ba U ^n(p, c(e)ro). 
This completes the proof of Theorem l5.ll □ 



6 Injectivity radius of null cones 

We turn our attention now to null cones within foliated Lorentzian manifolds. 
Our main result (Theorem l6 . 1 I below) provides a lower bound for the null injec- 
tivity radius under the main assumption that the exponential map is defined 
in some ball and the null conjugate radius is already controled. Hence, con- 
trary to the presentation in Section [3] our main assumption (see (A3') below) is 
not directly stated as a curvature bound. However, under additional assump- 
tions, it is known that the conjugate radius estimate can be deduced from an If 
curvature bound, so that our result is entirely relevant for the applications. 

Indeed, in a series of fundamental papers fl7l [TBI [19| , Klainerman and 
Rodnianski assumed on an L 2 curvature bound and estimated the null conjugate 
and injectivity radii for Ricci-flat Lorentzian (3 + l)-manifolds. Our result in 
the present section is a continuation of the recent work [ 19 [ and covers a general 
class of Lorentzian manifolds with arbitrary dimension, while our proof is local 
and geometric and so conceptually simple. 

We use the terminology and notation introduced in Sectional In particular, 
a point p e M and a reference vector field X are given, and N~ denotes the past 
null cone in the tangent cone at p. The null exponential map exp^ : B^(0, r) — » M 
is defined over a subset of this cone, 

B£(0,r):= Bt(0, r) n N~, 

and allows us to introduce the (past) null injectivity radius Null Inj (M,p, T). 
We also set 

03~(p,r) := exp^(0,r)). 

We consider a domain Q C M containing some point p on a final slice Eo 
and foliated as 

Q= [J E f , peE . (6.1) 

tsi-i,o] 
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We assume that there exist positive constants Kq, K\, K2 such that 

e~ K ° <n<e K ° in Q, (Al) 

\&Tgk < Ki in Q, (A2) 

the null conjugate radius at p is ro (at least) and the null exponential map 
satisfies 

e~ Kl gT, P lu?(o,ro)^ (e x Pp)*(gT ls«(0,ro)) ^ *** Ib?(Pa) ( A3 ') 
and, finally, there exists a coordinate system on the initial slice £_i such that the 
metric g is comparable to the n-dimensional Euclidian metric g# in these 
coordinates: 

e" Ko < g |e_!< e Ko gE' in B E _ 1/F (p, r ). (A4') 
We refer to K2 as the effective conjugate radius constant. 

Theorem 6.1 (Injectivity radius of null cones). Let Mbe a differentiable (n + 1)- 
manifold, endowed with a Lorentzian metric g satisfying the regularity assumptions 
(Al), (A2), (A3'), and (A4') at some point p and for some foliation ft.ty . Then, there 
exists a positive constant z'o depending only upon the foliation bounds Kq,Ki, the null 
conjugate radius ro, the effective conjugate radius constant K2, and the dimension n 
such that the null injectivity radius of the metric gatp satisfies 

Null Inj g (M,p,T)>z . 

It is interesting to compare the assumptions above with the ones in Section[3] 
Assumptions (Al) and (A2) are concerned with the property of the foliation and 
were already required in Section[3] 

Assumption (A3') should be viewed as a weaker version of the L°° curvature 
condition (A3). Recall that, under the assumptions of Theorem 13.11 which 
included a curvature bound, an analogue of (A3') valid in the whole of O was 
already established in l|4.8t . It is expected that (A3') is still valid when the 
curvature in every spacelike slice is solely bounded in some W space. 

Indeed, at least when the spatial dimension is n = 3 and the manifold is 
Ricci-flat, according to Klainerman and Rodnianski lUTl lTBj Assumption (A3') 
is a consequence of the following L 2 curvature bound 

HRmjflbp^) < K' 2 (6.2) 

for some constant K' 2 > 0. 

Assumption (A4') concerns the metric on the initial hypersurface and is 
only slightly stronger than the volume bound (A4). Furthermore, according to 
Anderson [1J and Petersen [22J the property (A4') is also a consequence of the 
curvature bound 

HRmyllL-pu^) < K' 2 (6.3) 
for m > n/2 and some constant K' 2 > and a volume lower bound at every scale 

r"" Vol^ (S Eo (p, r)) >v , re (0, r ]. (6.4) 

In summary, by combining Theorem 16.11 above with the results in IIT91 H6l 
we conclude: 
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Corollary 6.2 (Einstein field equations of general relativity). Let (M,g) be a 
Lorentzian (3 + \)-manifold satisfying the vacuum Einstein equation 

Kic g = 0. (6.5) 

Suppose that near some point p e M there exists a foliation Q of the form \6.1\ satisfying 
Assumptions {A\)-{A2) and such that the 1? curvature assumption (|6.2b holds on the 
initial spacelike hypersurface Z_i. Then, there exists a positive constant z'o depending 
only upon the foliation bounds Kq,K\ and the curvature bound K' 2 such that the null 
injectivity radius satisfies 

NxxllIn} g (M,p,T)>i . 

Proof of Theorem [67TI Step 1. Localization of the past null cone N~(p) between 
two flat null cones. Assumption (A3') provides us with a bound on the null 
conjugate radius, we need to control the null cut locus radius. We proceed as 
in Section 3] and introduce coordinates near the point p such that x a (p) = 0. 
Precisely relying on Assumptions (Al), (A2), and (A4'), we determine the 
coordinates X = (x a )sothatx° = t and the spatial coordinates {x>) are transported 
(via the gradient of the function t ) from the coordinates prescribed on the initial 
slice ZLj. Then, the Lorentzian metric reads g = -n 2 dt 2 + gjj dx'dxi and satisfies 
for some Co, C\ > 

t5- < n 2 < C , -J- 6, 7 < gij < Ci 6i jr (6.6) 

for all —ro < t < and (x 1 ) 2 + . . . + (x") 2 < (ro) 2 , and in these coordinates 
the reference Riemannian metric gj is comparable to the (n + l)-dimensional 
Euclidian metric gE := dt 2 + (dx 1 ) 2 + . . . (dx n ) 2 : 

1 

— g E < gr < Ci g E . (6.7) 

Denote by 23 %{q, r) the Euclidean ball with center q and radius r. Note that these 
inequalities holds within a neighborhood of p in O. The forthcoming bounds 
will hold in a neighborhood of the past null cone only. To simplify the notation 
we set 

1 1 
C ° " Co' Cl := C? 
In each time slice of parameter value t - a we introduce the n-dimensional 
Euclidian ball with radius b 

A a <b := [t = a, (x 1 ) 2 + ... + (x") 2 < b 2 } c E„, 

which is centered around the point p' with coordinates (a, 0, ••• ,0). We also 
define A" >b , J 5 ^^,. . . similarly. 

For any point q in a slice Z to satisfying -r < t < and x 1 ^) 2 H h x"(^) 2 < 

c 2 t 2 we consider the line (for the Euclidian metric) connecting q top: 

yW = (z,j-x\q),--- ,-*"(?)), TG[i ,0]. 
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This is a timelike curve for the Lorentzian metric g, since 



|y'(x)| 2 = -n 2 + gij ^^<- Co + Cl <0, 
to ro 

which shows that 

^ncrfp), te(-r ,0). 
On the other hand, we claim that the larger Euclidian cone A 1 „ ... contains 

° <Ci |t| 

the null cone, in other words 

4 CiM c (*T(p) U r(p)) c , f 6 (-d r ,0). 

Indeed, arguing by contradiction we suppose there exist a time to £ (— C\ To, 0) 
and a point q e -A> c ( connected to p by a causal curve y = y(s) with y(0) = p. 
After reparametrizing (in time) the curve is necessary we can assume that 
y(x) = (z,x ! (z)) for some ? < x < 0, as long as the point y(x) lies in the 
coordinate system under consideration. For this part of the curve at least we 
have 

, ,,o o dx 1 dx> 

o>l/l 2 = - 2 + ^^ 

which by d6.6b implies that f^-) 2 +. . -+{^&) 2 < C\ Cq, Therefore, after integration 
we find 

(x 1 ^) 2 + ■ ■ ■ + x n (q) 2 ) V2 (t' ) < VQ Ci t' < VQd r < r . 

Hence, we can choose t' Q - to, the whole curve lies in the system of coordinates, 
and is parametrized in the form y(x) = (x, x(x)), (x e [to, 0]). Moreover, we have 
|x(f )| < VCi Co | to I < Ci |fo|, which contradicts our assumption q e A^ f . 

In conclusion, we have localized the slices of the past null cone within 
"annulus" regions: 

7T(p) n L t c A\ Ci %cm , t e [-ci r , 0]. 



Sfep 2. The past null cone N~(p) viewed as a graph with bounded slope. We 
now obtain a Lipschitz continuous parametrization of the null cone. For any 
fixed q e A~^° we consider the vertical curve passing through q: 

y q (x) = (x, x 1 (q), ■■■ , x n {q)), x e [-cir , 0]. 

By Step 1 we know that there exists t 9 such that y«(T 9 ) e N~(p). Moreover, is 
unique since N~(p) is achronal, and this defines a map 

such that F(^) = y q {T q ). It is obvious F(-cir , 0) = p. 
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We claim that the map F is Lipschitz continuous with Lipschitz constant less 
than C\, as computed with the Euclidean metric E. Namely, by contradiction, 
suppose that |F(^i)-F(^ 2 )Ie > Q \q1-q2\E for some qi,qi g .A~T° , then by H6.7H in 

Step 1, F((Ji) would be chronologically related to F(q 2 ) and this would contradict 
the fact that N~(p) is achronal. Moreover, from Step 1 it follows that 

F(^)3N-(p)na3 E (p,c3r ). 

Step 3. Constructing an homotopy of curves on the null cone. Suppose that 
J\, 72 are two (past) null geodesies from p satisfying 

yi(0) = y 2 (0), |yi(0)| T = l/ 2 (0)| r = 1, 
71 (si) = 72(82). 

We claim that max(si,s 2 ) > c^tq, which will establish the desired injectivity 
bound by setting z'o = c^tq. 

We argue by contradiction and assume that max(si,s 2 ) < c^Tq. Taking into 
account Assumption (A2) and applying exactly the same arguments as in Step 1 
of Section S] we see that the gr-lengths of the curves y\, y% satisfy 

L( yj ,g T )< Sj e cc ^ <cl +3/i r (j = 1,2). 

By Step 1 of the present proof we know that the Euclidean lengths of y\, y 2 
satisfy 

L( 7j ,g E )<cl +1/i r (/ = 1,2). 

In particular, y\, y 2 c N~(p) n 23e(p, c^yq) and we can thus concatenate the curve 
71,72 an d obtain 

y ■- y~ 2 l U 71 : [0,sj + s 2 ] -> W{p) n B £ (p,c^ ). 

Since F(A~^°) 3 N~(p) n I>E(p,Cjro), there exists a smooth family of curves 
o t : [0,S! + s 2 ] — > N~(p) such that 

01 = 7, CJ = p, 

CT £ (0) = (T £ (Sl + S 2 ) = p, £ £ [0, 1]. 

Specifically, we choose 

o £ (s) := F(eF"V(s)), 

where the multiplication by £ is defined by relying on the linear structure of 
•A~ c Ji r ~ Br"(0, c\ro). Equivalently, by setting x'(s) = x ! (y(s)) we have the explicit 

formula 

o E (s) =f(~ c\r , ex 1 (s), ■■■ , £x"(s)). 
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It is clear that the Euclidean and gr-lengths of o £ satisfy 

L(o £ ,g E ) < e(l + C 1 )L(y,g £ ) < c\ +ll% r , 
L(a e ,g T )<4 +5/8 r . 

By Assumption (A3') on the null conjugate radius, we can lift to the null 
cone of the tangent space T p M the continuous family of loops a £/ and we obtain 
a continuous family of curves a £ defined on [0, Si + S2] such that 

ct £ (0) = 0, L(a Er gT,p) < c\r Q . 

Observe that the property L(o E , g T , p ) < c^r < r guarantees the existence of this 
continuous lift. By continuity all of the curves a £ are loops containing 0. As 
observed earlier in the proof for the case of bounded curvature, o\ consists of 
two distinct segments which clearly can not form a closed loop and we have 
reached a contradiction. □ 



7 Injectivity radius of an observer in a Lorentzian 
manifold 

Main result 

We are now a in a position to discuss and prove Theorem 11.11 stated in the 
introduction. As we have seen in the proof of the previous section, once the 
injectivity radius is controled, one can construct a foliation satisfying certain 
"good" properties. On the other hand, the concept of injectivity radius is clearly 
independent of any prescribed foliation. As this is more natural, we will now 
present a general result which avoids to assume a priori the existence of a 
foliation. This will be achieved by relying on purely geometric and intrinsic 
quantities and constructing coordinates adapted to the geometry. Such a result 
is conceptually very important in the applications. The result and proof in this 
section should be viewed as a Lorentzian generalization of Cheeger, Gromov, 
and Taylor's technique [10|, originally developed for Riemannian manifolds. 

Let (M, g) be a differentiable (n + l)-manifold endowed with a Lorentzian 
metric tensor g, and consider a point p e M and a vector T e T p M with 
gp(T,T) = -1. That is, we now fix a single observer located at the point p. 
As explained in Section [2] the vector T induces an inner product gj — {, )j on 
the tangent space T p M. We assume that the exponential map exp^ is defined 
in some ball Br(0, tq) c T p M determined by this inner product, which is of 
course always true in a sufficiently small ball. Controling the geometry at the 
point p precisely amounts to estimating the size of this radius ro where the 
exponential map is defined and has some good property. We restrict attention 
to the geodesic ball ®r(p, ro) := exp f) (Br(0, ro)); recall that these sets depend 
upon the vector T given at p. 
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As explained in the introduction, by g-parallel translating the vector T at 
p along a geodesic y from p, we can define get a future-oriented unit time- 
like vector field T — y defined along this geodesic. To this vector field and 
the Lorentzian metric g we can associate a reference Riemannian metric gr 
along the geodesic. In turn, this allows us to compute the norm |RmJr of the 
Riemann curvature tensor along the geodesic. 

Of course, whenever two such geodesies y, y' meet away from p, the corre- 
sponding vectors T y and Ty are generally distinct If we consider the family of 
all such geodesies we therefore obtain a (generally) multi-valued vector field 
defined in the geodesic ball Ht(p, Tq). We use the same letter T to denote this 
vector field. In turn we can still compute the Riemann curvature norm |Rm ? |x 
by taking into account every value of T. 

The key objective of the present section is the study of the geometry of the 
local covering exp p : Bj(0, Yq) — » 23j(p, yq) and to compare the Lorentzian metric 
g defined on the manifold M with the reference Riemannian metrics gj. As we 
will see in the proof below, it will be convenient to pull the metric "upstairs" on 
the tangent space at p, using the exponential map. Indeed, this will be possible 
once we will have estimated the conjugate radius (in Step 1 of the proof below) 
and will know that the exponential map is non-degenerate on Bj(0, Yq). Pulling 
back the Lorentzian metric g on M by the exponential map we get a Lorentzian 
metric g = exp*g defined in the tangent space, on the ball Bt(0, yo). We use the 
same letter g to denote this metric. Then, the geometry in the tangent space 
is particularly simple, since the g-geodesics on M passing through p are radial 
straightline in Bt(0, Yq). 

A third view point could be adopted by restricting attention within the cut- 
locus from the point p, and by imposing the curvature assumption within the 
cut-locus only 

We are in a position to prove the main result of the present paper that was 
stated in Theorem ll.il 

Pyoo/ of TheoYem [TT1 After scaling we may assume that Yo = 1, and so we need 
to show 

Inj g (M,p, T) > c(n) Vol,(25 r (p, c{n))). (7.1) 

Step 1. Estimates for the metric gj and its covariant derivative. Let Eq = T, 
Ei, ■■ ■ ,E n be an orthonormal frame at the origin in T p M for the Lorentzian 
metric g. By g-parallel transporting this basis along along a radial geodesic 
y = y(r), satisfying y(0) = 0, |y'(0)|r = 1, we get an orthonormal frame defined 
along the geodesic. We use the same letters E a to denote these vector fields. 
Since ^ 

— (E a ,Efs) g = 0, 

we infer that 

\Ej\j = |Ej|g = 1 along the geodesic. 
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The same argument also implies 

l/Mlf = ly'(0)| 2 T = 1, \y'(r)\ 2 g = \V(0)\ 2 g = 1, (7.2) 

and y'(r) = c a E a (r) with constant (in r) scalars c" and £, \c a \ 2 = |y'(0)|r = 1. We 
used here that, by definition, y' is ^-parallel transported. 

Let V - a a (r)E a {r) be a Jacobi field along a radial geodesic y = y(r), with 
V(0) = and \V (0)|r = 1. Then, the Jacobi equation takes the form 

(a a )"(r) = -(E a ,R(E p ,E r )E 6 ) T cPc 5 a?(r). 

Since 

-2 £ («; 2 + 4) < < 2 + ««) * 2 E (< 2 + 4 

we obtain |V'(r)| T < e r and thus \V(r)\ T < (e r - 1). 

By substituting this result into the above formulas, the estimate can be im- 
proved again. Indeed, by computing and estimating the second-order deriva- 
tive £ Yj a a' a a a as we did for the Jacobi field estimate of Section|4j we can check 
that 

r-C{n)r 2 <(J^\a a \ 2 {r)) ' < (e r - 1) along the geodesic. 

Denote by go, gj,o the Lorentzian and the Riemannian metrics at the origin 
(which are nothing but the metrics at the point p), and let y°,...,y" be Cartesian 
coordinates on Br(0, 1), with , ^) ga (0) = t] a p. Assuming that the radius 
under consideration is sufficiently small so that (1 - C(n) \y\) < 1 we conclude 
from the Jacobi field estimate that the exponential map is non-degenerate and 
that the metric along the geodesic are comparable. In turn, since this is true for 
every radial geodesic, we can define the pull back of the metric to the tangent 
space and the conclusion hold in the whole ball Bj(0, 1), that is 

(l-C(n)\y\)gT,o<gT,y<(l + C(n)\y\)g T , , yeB r (0,l). (7.3) 

By construction of the metric gj we have V ?r — V g = V g T * T (schematically) 
and VT(0) = 0, and it is useful to control the covariant derivative too. To this 
end, write the radial vector field as 

with = 1 (as stated already in O). Using that \VT\ 2 = V a T s V g rM gf 
and computing the derivative of |VT|| along radial geodesies, we find 

^-\VT\ 2 T < C(n) \VT\l + 2<Vi VT,VT)t. 
dr * 

By using that 

V7 nn „ r d 9 . Id if d 

VsT = 0, — ,— ]= — + ^r—, 

* dr dy a r dy a r 2 dr 
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we obtain 

Va V a V = — V a V + R(— , — )V , 

Tr w r ? dr dy a 

and therefore, thanks to the curvature assumption, 



d ,. 2. w ., 

"r- 



f- |VT|2 < — |VT|| + C(n) |VT| 3 r + C(n) |VT| 



This implies the following bound for the covariant derivative 

|VT| r (y) < C(n) \y\, |y| < 1/C(n), (7.4) 
which also provides a bound for the difference V ?T - V ? . 

Step 2. Estimate of the injectivity radius of g on £>t(0, c{ri)). 

Since the curvature on Br(0,l) is bounded and that | V ?r -Vg|?. < C(n) = l/c(n) 
on the ball £>r(0, c(n)) we can follow the argument in Section|4jand bound from 
below the conjugate radius for any point in the ball Bt(0, 3c(n)/4). 

Next, given any point y 6 Bj{Q, c(n)/2), let y\ and j2 be two geodesies which 
meet at their end points and have "short" length with respect to the metric gj 
(or gT,o)- By using the linear structure on Bj(0, 1) (a subset of the vector space 
TpM) we can construct an homotopy of the loop j\ U y~ l to the origin, such that 
each curve have also "short" length for the metric gj. By lifting the homotopy 
to the tangent space T y Bj(0, 1) and by relying on the conjugate radius bound, 
we reach a contradiction as was done in Section[4] 

In summary, there exists a universal constant C(«) = l/c(«) depending only 
on the dimension such that the injectivity radius at each point y of Bj(0, c(n)) is 
bounded from below by 4c(n). Moreover, by the Jacobian field estimate again, 
we can prove the ball B T , P (0, c(n)) c T p M defined by the Euclidean metric gr, P is 
covered by exp v (Br /y (0, 3c(n))), where B T , V (0, 3c(n)) c T y T p M is a ball of radius 
3c(n) defined by metric gr,y, and any two points in Bt jP (0, c(n)) can be connected 
by a g geodesic which is totally contained in 6^(0, 2c(n)). Further arguments 
are now required to arrive at the desired bound i7.1) . 



Step 3. Riemannian metric induced on Bj(0,2c(n)). Consider a geodesic y 
satisfying y(0) = and y'(0) = -T, and let us set 

y(c(n)/2) =: q, z := d g (; q) - d g (q, 0). 

Then, by following exactly the same arguments as in the main proof of Section|5l 
we construct a normal coordinate system (of definite size) such that g = -dx 1 + 
gij dx'dx 1 and gw = dr 2 + gijdx l dx ! , and such that the corresponding reference 
Riemannian metric satisfies the following properties: 

(i) (1 - c{n)) g N <g T <(l + c(n)) g N , 

(ii) gN has bounded curvature (< C(«)) (see J5.5|l ), and 
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(iii) for any fixed yo 6 B T (0, c(rij) the distance function d gN (y , ) 2 is strictly 
g-convex on the ball Bt(0, 2c(n)) and, more precisely, 

(2 + c(n)) g N > V 2 dJ(y , •) > (2 - c(n)) g N on B r (0, 2c(n)) 

for any y e B T (0, c(n)). 

The Hessian of the distance function (defined by the Riemannian metric giv) is 
naturally computed using the covariant derivatives defined by the Lorentzian 
metric g. 

Step 4. Suppose that pi, •■ ■ , p;v are distinct pre-images of pin the ball £>r(0, c(tt)). 
We claim that any p' e 23r(p, c(n)) has at least N distinct pre-images in Br(0, 1), 
and refer to this property as a "lower semi-continuity" property. 

Generalizing the terminology in 1 10 1, we use the notation a ~ b when two 

(gTflA) 

curves a, b defined on M and with the same endpoints are homotopic through 
a family of curves whose lift have gxo'lengms < A. Relying on the lift and the 
linear structure, we see that, for any curve £, starting from p with (after lifting 
through 0) gr^-length A < 1, there exists a unique g-geodesic y% (with the 
same end points as defined on M such that £ ~ y^. This fact establishes a 

(gTflA) 

one-to-one correspondence between the following three concepts: 

(i) equivalence class of curves through p with g^crlengths < 3c(n), 

(ii) radial geodesic segments of gr^-lengths < 3c(n), and 

(iii) points in the ball B T (0, 3c(n)) c T p M. 

Let a be a g-geodesic connecting p to p' in Sr(p, c(n)). Observe that the 
images of the lines Op, by the exponential map, a, = exp p (Op ( ), are distinct 
geodesic loops through p. Denote by cr, the lift of a; U a through 0, and denote 
by p' the end point of <j,. Then, it is clear that all the points p' t (i = 1, • • • ,N) are 
the pre-images of p' in Bt(0, 1/2). We claim that they are distinct. 

Indeed, assuming that p' = p' for some i + i, we would find crUcr; ~ ffUcr,-, 

> gTflM"-) 

which gives 

Oi ~ a -1 U a U Oj ~ a~ l UaUff; ~ o;. 

g T p,3c(n) gT,a,3c(n) g Tfl ,3c(n) 

This would imply a, ~ a, and, therefore, p, = p ; , which is a contradiction. In 

gT,oMn) 

short, this argument shows that the "cancellation law" holds for the homotopy 
class of "not too long" curves. 

Step 5. Suppose that there exist two distinct g-geodesics y\ : [0, h] — > M and 
yi '■ [0, h] —> M satisfying 

yi(0) = y 2 (0) = p, |y'(0)| 2 r = |y'(0)| 2 r = 1, 

and meeting at their endpoints, that is: yi(h) = 72(b)- Then, let I :=l\+ I2 and 
y :- 1 U yi : [0/ ~~ * M. Our aim is to prove that 

l>c(n)Vol g ('BT(p,c(n))), 
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which will give us the desired injectivity radius. 

From the loop y we define a map n y : B T (0, c{n)) — > B T (0, 2c(n)) as follows: 
for any y e B T (0,c(n)), the point n y {y) is the end point of the lift exp p (Oy) U y 
(through the origin). If one would have n y (y) = y then by the cancellation law 
established in Step 4, we would have y ~ 0, which is a contradiction. So, 

g T ,o,2c(n) 

the map n y has no fixed point. 

Without loss of generality we assume that I < c(nf. Let N = [c(n) 3 /l] be the 
largest integer less than c(n) 3 /l, and let us use the notation 2y = y o y, etc. 

Claim. The classes [y], [2y], • • • , [Ny] are distinct homotopy classes for the 
relation ~ . 

g T ,oAn) 2 

If this were not true, then by the cancellation law we would have [jy] ~ 

gTfiA") 2 

for some 1 < j <N. We already know that all n l y is defined from B T (Q, c(n) 2 ) to 
Bt(0, c(n)) for i < j. Since for any y e Bj(0, c(n) 2 ) we have 

exp (Oy)Ujy ~ exp (Oy), 

r gTfiAn) v 

which implies that n y = id. We use here the notation n y := n y o n y , etc. 
Then, we define a function u : B T (0, c(n)) — > ]R by 

u(y) = d|(0, y) + d|(0, n 7 y) + ■■■ + d|(0, Y y). 

Since nj, = id, it is easy to see u(n y y) = u(y) for any y e B T (0, c(n)). That is to say, 
u is 7ij,-invariant. By Step 3, u is strictly g-geodesically convex on B T (0,c(n)). 
More precisely, since for any g-geodesic £ : [0,s ] — > B r (0, c(n)), 71^4 are still 
^-geodesies in B r (0, c(n)), and 

^u(m = V 2 d|(0, -)(5'(s), «'(s)) + • • • + V 2 d|(0, ■)(dn! y -\ {s) (e(s)Un y ~\ (s) (e(s))) 

>g(e(s),e( S ))>o. 

Observe that 

2 2/c(n) 3 , ;'c(n) 2 
u lB r (o,c(«))' > ;(1 - c(n)) (c(n) — ) > — — , 

and 2 
u{Q)<j{jlf<jc{nf< ] -^-, 

so the minimum of function u over Br(0, c{n)) is only achieved at at an interior 
point, say yo e B T (0,c(n)). Then by n y invariance of u, we have u(n y yo) = 
u(y ) < jc(n) 2 /2, and this implies n y (yo) e B r (0,c(n)). By the injectivity radius 
estimate at yo 6 (T p M,g), there exists a g-geodesic connecting yo to n y (yo), 
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which is contained in Bj iV (0,2c(n)). By using the strong g-geodesic convexity 
of u, we conclude that n y yo = i/o- This contradicts the fact that n y has no fixed 
point, and the claim is proved. 

Step 6. The pull back of the volume element of g is the same as the one of gj. 
By combining this observation with our results in Steps 4 and 5 we find 

Vol, T (B r (0,l)) > C -^fvol g CB T (p,c(n))), 



which gives 

Vol g (3 r (p,c("))) 
VoV(B T (0,l)) 

The proof of Theorem 1 1.1 1 is completed. 



l> ~ C{n) VoMBHO,!)) Mn)V 0W M»)))- 



8 Volume comparison for future or past cones 

In Riemannian geometry, under a Ricci curvature lower bound, Bishop-Gromov's 
volume comparison theorem allows one to compare the volume of small and 
large balls in a sharp and qualitative manner. Let us return to Step 2 of Sec- 
tion |5j where we introduced the index form associated with the synchronous 
coordinate system on time-like geodesies. By noticing that the index form is 
symmetric and that Jacobi fields minimize the index form (in some sense), we 
can extend the method of proof of the index comparison theorem. However, 
in a general Lorentzian manifold, since the index form we needed (without 
imposing a restriction on the geodesies) is non-symmetric, we need to adapt 
the method of the index comparison theorem, as follows. 

Theorem 8.1 (Volume comparison theorem for cones). Let (M,g) be a globally 
hyperbolic, Lorentzian (n + l)-manifold. Fix p e M and a vector T e T p M with 
g p (T,T) = — 1, and suppose that the exponential map exp p is defined on the ball 
Bt(0, r ) c T p M (determined by the reference inner product g T at p). Suppose also that 
the Ricci curvature satisfies on Urip, r o) 

Ric g (V, V) > -n K 2 1 1 V||| for all time-like vector fields V. 

Then for any < r < s < ro the inequality 

Vol g ( 9C(p,r)) > Vol K2 (B(r)) 
VolgPCfaa))) ~ Vol K2 (B(s))) 

holds, with "JC(p, r) := exp p (FC(p, r)) and 

FC(p,r) := (0 < \V\ gTfi < r , \V\ 2 g < 0, (T, V) gTfl < oj 



33 



and Voljc 2 (B(r)) is the volume of the ball with radius r, B(r) (analogous to Urip, r ) c M), 
in the simply-connected Lorentzian (n + l)-manifold with constant curvature K 2 (that 

IS, Ra(iyt> = ~K 2 (gaygp - gadgfty))- 

More generally, if E is a subset in unit sphere S" such that \V\ g < Q,(T, V) g < 
for all V e E, then the inequality 

Vol g (?C z (p,r)) > Vol K2 (B(r)) 



Vol g (?C z (p,s))) ~ Vol K2 (B(s))) 
holds with 3Cz(p, r) := exp p (FCz(p, r)) and 

FC z (p,r):=[VeFC(p,r)/^-eL}. 

This result will be used shortly to control the injectivity radius of null cones, 
but is also of independent interest. For definiteness we state the result for future 
cones. 

Proof. Let y : [0,s ] — > M be a future-oriented time-like geodesic satisfying 
y(0) = p and |y'(0)| ?T = — 1. We are going to use the standard technique to 
compute the rate of change of the volume element along y. Given si e (0, so) 
assume that any point in (0, Si] is neither a conjugate point nor a cut point with 
respect to p. Let v = y'(s\), V\, v 2 , • • • ,v n be an orthonormal basis at y(si) with 
respect to g y ( Sl )- Let also V a be the Jacobi field defined on [0,si] and satisfying 
V a (0) = and V a (s{) = v a . Clearly V = (s/sj) y' , and the vectors Vi and VyVj 
are orthogonal to y' for all i > 1. 

Consider the Jacobian of the exponential map <p(s) := /(dexp y(s) ), which is 
given by 

2 _ \y'(s)AV,(s)A---AV n (s)\ 2 g 
<P(S) " s 2 "\y'(0)AV[(0)A---V' n (0)\ 2 g 
Denote by (pK 2 (s) the corresponding quantity in the simply connected Lorentzian 



(n + l)-manifold with constant curvature -K 2 . Define the index form 



I S1 (X,Y) := j 1 ({V r X,V r Y) g - Rm g (y',X,y',Y))ds, 



where X, Y are vector fields along y and Rm ? (y', X,y',Y) := -{Rm g (y', X)y' , Y) g . 
Observe that I Sl is symmetric in X, Y. It is easy to see 



d_ 

ds 



log<p 2 = Y(V' i (s 1 ),V i (s 1 )) g -^ 

-Si ' . ■ S\ 

i 

E2n 



Let E,(s) be the parallel transport of i>; along y. Since there are no conjugate 
points along y, the Jacobi field minimizes the index form among all vector fields 
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with fixed boundary values. This is the same as in Riemannian geometry. The 
reason is that the length of time-like geodesic without conjugate points is locally 
maximizing among all nearby time-like curves with the same end points. Let 

Vi(s) = ffiEKs), then I Sl (V t , V,-) < l Sl {V h V t ) and 

£u log < - £ f ¥^<W/. Eu y%) - K 2) 
as ° ((pK 2 ) z *-f Jo (smhsi) 2 

= _ p jsinh_s)V _ ^ ^ 

Jo (sinhsi) 2 

The following is a simple but very important observation due to Gromov, 
which we now extend to a globally hyperbolic Lorentzian manifold. Let A 
be the star-shaped domain (with respect to 0) in T p M, such that exp p : A n 
f>r(0, J"o) is a diffeomorphism on its image and the image of dA D Bj(0, ro) is 
set of cut locus (in "Bjip, ro)). Let xa be the characteristic function of A. Since 
<p(s)/<pK 2 (s) is decreasing in s we see that Xa<p/<Pk 2 is also decreasing in s. Now, 
we get two functions on Bj(0, ro), whose quotient is decreasing along radial 
geodesies. Observe that M is globally hyperbolic, so any point in 3 r C(p, ro) 
is connected to p by a maximizing time-like geodesic. This also implies that 
the integration of xa<P over Bj(0,s) gives the volume Vol^(!Br(P/S)). Then, 
by integrating xa<P and <Pk 2 over £>r(0, s) and after a simple calculation we 
deduce that Vol g (3 r C(p,s))/VolK 2 (B(s)) is decreasing in s. The case of the ratio 
Volg(9 r Cj:(p, s))/Voljc 2 (B(s)) is similar. The proof of the theorem is completed. □ 

We are now in a position to prove : 

Corollary 8.2 (Injectivity radius based on the volume of a future cone). Let M 
be a manifold satisfying the assumptions in Theorem \T7l\ and assumed to be globally 
hyperbolic, and let T e T p M be a reference vector. Let Lbea subset in the unit sphere 
S n included in the future cone N*. lf\o\ g (3 : C L z{p, Tq)) > v Q > 0, then the inequality 

^) S W,P,T) Vo 

Vo ^ C(I V 

holds, where JC^p, r ) := exp f) (FC p (r )) with 

FC p (r ) := jo < |V| r < r , (T, V) T < 0, \V\ 2 g < 0, 6 zj, 

and the constant c(E) depends only on the distance (measured by T) of E to the null 
cone. 

Proof. First, we recall there is a constant C(E) depending only on the distance of 
E to the null cone, such that Ric()/, y') > -C(E)|y'|^ r for any time-like geodesic y 
with y'(0) e E. By the volume comparison theorem for future cone established 
in Theorem l8.1l we have 

Vol g (?C z (p,c(n)r )) 



Vol g (?C z (p,r )) 



> C(E), 
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and combining this result with Theorem ll.il the corollary follows. □ 

9 Final remarks 
Regularity of Lorentzian metrics 

Following the strategy proposed in the present paper, we now "transfer" to the 
Lorentzian metric the regularity available on a reference Riemannian metric. 
Clearly the regularity obtained in this manner depends on the way the refer- 
ence Riemannian metric is constructed. The interest of our approach below is 
to provide a simple derivation: using harmonic-like coordinates for the Rie- 
mannian metric we see immediately that the Lorentzian metric has uniformly 
bounded first-order derivatives. For the optimal regularity achievable with 
Lorentzian metrics we refer to Anderson (3). 

Proposition 9.1 (Regularity in harmonic-like coordinates). Under the assump- 
tions and notation ofTheorem \l.l\ define 

Vol g CB T (p,c(n)r )) 

f i := C W r °' 

f o 

where c(n) is the constant determined in this theorem. Then for any e > there exist 
a constant C\{n,e) with lim t ^ c i( n , £') - and a coordinate system (x a ) satisfying 
x a (p) = and defined for all (x ) 2 + (x 1 ) 2 + . . . + (x") 2 < (1 - e) 2 r\, such that 

\g a? -r\ a p\<ci{n,E), 
n\dgap\<ci(n,e), 

where rj a p is the Minkozvski metric in these coordinates. 

Proof. By scaling we may assume r\ = 1. By Step 1 in the proof of Theorem ll.il 
we know that the Riemannian metric gj is equivalent to the Riemannian metric 
gr,o on £>r(0,4c(«)). By considering a lift and using again the results in Step 1 
this implies 

3 T (p, c(n)) c S T ( 9/ 3c(n)) q 6 2 T (p, c(n)). 

Applying the same argument as in Theorem ll.il we deduce that the injectivity 
radius of any point in T>t(p, c{n)) is bounded from below by c(n). As in Step 3 
in the proof of Theorem ll.il (or in Step 2 of Section [5), we see that there exists 
a synchronous coordinate system (y a ) - (t, y') of definite size around p such 
that the metrics g = -dz 2 + gij dy'dyi and gN = d% 2 + gij dy'dyi (the Riemannian 
metric constructed therein) satisfy the following properties on the geodesic ball 
•BT(p,c(n)): 

(a) (1 - c(n)) g N < gr < (1 + c(n)) g N , 

(b) gN has bounded curvature (< l/c(n)), 
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(c) |t| + fa + |V 2 t| n < l/c(n). 

(In particular, this implies |V ?Ni g|jv < l/c(n).) Since the volume Vol gN (S r (p,c(n))) 
is bounded from below, it follows from |10| that the injectivity radius of gjv at 
p is bounded from below by c(n). By the theorem in |16| on the existence of 
harmonic coordinates, for any small e > there exists an harmonic coordinate 
system (x a ) with respect to the Riemannian metric gN such that £ a \x a \ 2 < (1-e) 2 
and for every < y < 1 

\gN A p ~ M < c i("' £ )> I^nI < l/c(n), I^nIo < l/c(n, £,y). 

In the construction of harmonic coordinates, we may also assume that |^- - 
JfUo <cx(n,£). 

Since |V ?N g|jv < l/c(n) and that, in these coordinates, |V^ N | < l/c(n), we have 
\dg\ < l/c(n). Finally, to estimate the metric we write \g a p - r\ a p\ v < c\(n, e) and 
\dg\ < l/c(n) and we conclude that \g a p - t] a p\ < + C\(n,e). The proof is 
completed. □ 

Pseudo-Riemannian manifolds 

Finally, we would like discuss pseudo-Riemannian manifolds (M, g) (also re- 
ferred to as semi-Riemannian manifolds). Consider a differentiable manifold M 
endowed with a symmetric, non-degenerate covariant 2-tensor g. We assume 
that the signature of g is (n\, n-i), that is, ti\ negative signs and «2 positive signs. 
Riemannian and Lorentzian manifolds are special cases of pseudo-Riemannian 
manifolds. Fix p e M and an orthonormal family T consisting of n\ vectors 
Ei/E2/'" ,E Hl e TpM such that (Ej,Ej) g = — 6«. Based on this family, we can 
define a reference inner product gj on T p M by generalizing our construction 
in the Lorentzian case, and by using this inner product we can then define the 
ball Bj(0, s) C T p M. By parallel translating E\, E2, ■■■ , E, h along radial geodesies 
from the origin in T p M, we obtain vector fields E\, E2, ■ ■ • , E ni defined in the tan- 
gent space (or multi-valued vector fields on the manifold). This also induces a 
(multi-valued) Riemannian metric gj as was explained before. 

The following corollary immediately follows by repeating the proof of Theo- 
rem ll.ll We note that the curvature covariant derivative bound imposed below 
is probably superfluous and could probably be removed by introducing a foli- 
ation based on certain synchronous-type coordinates, as we did in Section|5]for 
Lorentzian manifolds. On the other hand, to the best of our knowledge this is 
the first injectivity radius estimate for pseudo-Riemannian manifolds. 

Corollary 9.2 (Injectivity radius of pseudo-Riemannian manifolds). Eet (M, g) be 

a differentiable pseudo-Riemannian n-manifold with signature (m, nj), and letpeM 
and T = (Ei,--- ,E ni ) be a family of vectors in T p M satisfying g(E„ Ey) = -6,y. 
Suppose that the exponential map exp^ is defined on Bj(0, Tq) c T p M and that 

|Rn%| r < rf, \VRm g \ T < rf on B T (0, r ). 
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Then, there exists a positive constant c(n) such that 

Inj (M, p,T) Vol ? CB T (p, c(n) r )) 

^ c(n) , 

ro r\ 

where 23r(p, r) = exp p (Br(0, r)) is the geodesic ball at p with radius r. 

Proof. Without loss of generality we assume ro = 1 . In local coordinate system 
y a ,let 

£, =: Ef -j-J, E itt = E^g^, i = 1, . . ., m, 

then gT A p = g«js + 2 TZ-i EiaEip- By the same computations as in the proof of 
Theorem 1 1.1 1 we obtain 

IVEflr < 4t/ 
c(k) 

Igr - gr,ol + \g- r\\ < c(n) on the ball B T (0, c{n)), 

where rj a a := +b a p (a minus sign for a < n\, and a plus sign for a > n{). In view 
of the computations in [ 13 1 (Theorem 4.11 and Corollary 4.12) we deduce that 

\dg\ < r/c(n), where r 2 = (y 1 ) 2 H + (y") 2 . Since d| ro (y , y) = |y - yol 2 , we have 

for any point yo e B T {0, c(n)) 

V lp d2 gT,o(yo> •) ^ = on the ball B r (0, c(n)). 

Since the metric gr,o plays the same role as gu (cf . the proof of Theorem lLlL all 
arguments can be carried out and this completes the proof of the corollary. □ 
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